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1 Introduction 

It was proved by Weyl [8] in 1916 that the sequence of values of ain? is uni- 
formly distributed modulo 1, for any fixed real irrational a. Indeed Weyl's 
result covered sequences an'^ for any fixed positive integer exponent d. How- 
ever Weyl's work leaves open a number of questions concerning the finer 
distribution of these sequences. It has been conjectured by Rudnick, Sarnak 
and Zaharescu [6] that the fractional parts of ari^ will have a Poisson distri- 
bution provided firstly that a is "Diophantine" , and secondly that if a/q is 
any convergent to a then the square- free part of q is Here one says 

that a is Diophantine if one has 



a 

a 



for every rational number a/q and any fixed e > 0. In particular every real 
irrational algebraic number is Diophantine. One would predict that there are 
Diophantine numbers a for which the sequence of convergents Pn/q-n contains 
infinitely many squares amongst the qn- If true, this would show that the 
second condition is independent of the first. Indeed one would expect to find 
such a with bounded partial quotients. 

The Poisson property can be phrased in terms of a sequence of correla- 
tion functions. We shall be concerned in the present paper with the pair 
correlation function. For a real sequence 9 = {On)]° considered modulo 1 we 
define the pair correlation function by 

R{N,X]9) := N-'^4f{m < n < N : 119^ - 0n\ \ < XiV'^}, 

and if 6'„ = an^ we write Ra{N, X) in place of R{N, X; 0). If the sequence 6 
follows a Poisson distribution then we will have 

hm R{N, X;e) = X for all X > 0. (2) 

A''-+oo 
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The statement ([2]) is in general weaker than the Poisson condition. How- 
ever we know rather httle even about Ra{N,X). For the pair correlation it 
appears that one does not need the "nearly square-free" condition for the 
numerators q of the convergents a/q of a. We therefore make the following 
conjecture. 

Conjecture 1 If a is Diophantine, then 

RaiN.X) = N-^j^{m <n<N: \\a{m'^ - n^)\ \ < XN-^} X (3) 
as N ^ oo, for each fixed X > 0. 

Related conjectures are already mentioned in the works of Rudnick and 
Sarnak [5J and of Rudnick, Sarnak and Zaharescu [6] . 

We remark that if \a — a/q\ < l/(4g^) infinitely often, then ([3]) is false for 
every X G (1/4, 1/2). To see this one may take N = q and note that the pairs 
m and n = q — m for m < q/2 will satisfy ||Q;(m^ — n^)|| < 1/(4A^), whence 
Ra{N,X) > 1/2 + o(l). Thus some condition on rational approximations to 
a will clearly be necessary. A similar remark occurs in the paper of Rudnick 
and Sarnak [3]. 

Rudnick and Sarnak |3] were able to show that "almost all" a, in the 
sense of Lebesgue measure, satisfy (E]), but they remark that they are not 
able to provide any explicit value of a which does so. An alternative proof of 
this result was given by Marklof and Strombergsson [Aj. Our first result gives 
a third way to establish the "almost-all" property, but more importantly it 
allows us to construct values of a for which holds. 

Theorem 1 The statement ([3]) holds for almost all real a. Moreover there 
is a dense set of constructible values of a for which ([3]) holds. 

The second claim of the theorem deserves further comment. What we will 
do is to provide an informal algorithm, which, for any closed interval / of 
positive length, provides a convergent sequence of rational numbers belonging 
to /, whose limit a satisfies It could be said Rudnick and Sarnak were 
hoping for a more explicit construction, akin to that for Liouville numbers, 
for example. However, from a logical point of view there is no difference 
between our construction and that of other more familiar real numbers. The 
reader might also feel happier if we had given an explicit example of an 
admissible a, by displaying its decimal expansion; but since our construction 
provides a dense set of values, that would be uninformative. We can safely 
assert that 

a = 3. 14159265358... 
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satisfies 1^, but this will not help the reader's intuition! 

We note at this point that our proof of Theorem [T] provides slightly more. 
Indeed there is a positive constant rj (we may take rj = 1/200) such that 
for almost all a, and in particular for those a which we construct, we have 
RaiN, X)=X + 0(Ar-^) uniformly for A^"'? < X < iV''. 

Rudnick and Sarnak [5j proved that holds for the sequence (an'^)f^ 
for almost all a, for every d > 2. However our approach appears to work 
only for d = 2. 

In proving Theorem [1] we shall show that ([3]) holds for all a satisfying 
three conditions, which are explained in detail in §3. The first of these is 
that 



a 

a 

Q 



for every approximation a/q to a. It is of interest that this requirement is not 
quite as strong as ([I]). The second condition is roughly that if a„/g„ are the 
continued fraction convergents to a, then g„ is "almost odd and square- free" . 
It seems conceivable that one could adapt the proof to avoid this condition. 
The third assumption on a is that a„ does not lie in a certain small "bad" 
set B{qn), if n is large enough. One would conjecture that the sets B{q) are 
empty for all sufficiently large q. Thus in this approach it is the sets B{q) 
which are the real stumbling block in any attack on Conjecture 1. 

A related approach to Conjecture 1 has been investigated by Truelsen [7]. 
This is based on a hypothesis concerning the average value of the function 

Tlj{n) := #{(a, 6) G : a, 6 < M, ab = n} 

in short arithmetic progressions. Such a hypothesis is related to the condi- 
tion giving our bad sets B{q). Truelsen proves that his hypothesis holds on 
average, in a suitable sense. Our Lemma 3 is in a similar vein, but the two 
results are not directly comparable. 

Our second result gives partial support to Conjecture 1, by describing the 
behaviour of Ra{N,X) as X grows. 

Theorem 2 Suppose that a eM. and k > 1 satisfy 

a 



a 



1 

> 



for all fractions a / q . Then 

R^iN,X) = X + OiX''/^) + OiK'ihgN)-') 
uniformly in all the parameters, for 1 < X < logA^. 
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This result applies in particular whenever a is Diophantine. It shows 
that, in the limit as — ^ oo, the function Ra{N,X) is approximately equal 
to X for large X. Moreover we have the correct order of magnitude 

X<, RaiN,X) <«X 

as soon as X 1- 

For a fixed X < 1 we are unable to prove even that Ra{N,X) ^ 1 in 
general. However the method used to establish Theorem [2] can be adapted 
to yield some non-trivial upper bounds, of the form Ra{N,X) <^ X^ with 
9 > 0. Here we require a to be Diophantine, and (\ogN)~^ < X < 1 for a 
suitably small constant 6 > 0. 

In discussing Theorem |2] it is natural to examine the case a = a/q, which 
leads to consideration of congruences a(m^ — n"^) = r(mod q) with r small. 
Thus it would be interesting to know about the number of solutions u,v < N 
of uv = c(mod q), for a fixed c. During the proof of Theorem [2] we will use 
a result of Linnik and Vinogradov [3] which shows that 

d{k) <^s<Piq)q'^NHogN 

k=c (mod q) 

uniformly for (c, q) = 1 and q < N'^^^~^\ for any fixed 5 G (0, 1). (Here d(k) 
is the divisor function.) However for our problem we expect that the factor 
log N can be removed, and we make the following conjecture. 

Conjecture 2 For any fixed 6 G (0, 1) we have 

v) eN^ : u,v < N, uv = c (mod q)} <5 (l){q)q'^N'^ 
uniformly for (c, q) = I and q < N'^^^~^'^ . 

Unfortunately it appears that the techniques used by Linnik and Vinogradov 
do not work for the above variant of their problem. It is no coincidence that 
Conjecture 2 can be reformulated using Truelsen's function T^{n). 

In order to put Theorem [2] into context it may be helpful to record what 
one can say about arbitrary sequences 6. For this purpose it will be more 
convenient to use a weighted pair correlation function 



Ro{N,X-9):=N-' E 1- 
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Theorem 3 Let 9 = {9n)T arbitrary real sequence. 

(i) We have ([2]) if and only if 

\im Ro{N,X]e) = l+X forallX>0. (4) 

(ii) We have Ro{N, X; 9) > max(l, X) for all X >0. 

(ill) We haveRoiN,X + Y;9) < Ro{N,X;9) + RoiN,Y;9) forallX,Y > 0. 

Notice in particular that in part (i) we make no assumption about uniformity 
with respect to X in either of the hmits involved. We remark that Part (ii) 
can be strengthened slightly with a little more work. If X = [X] + ^, where 
[X] is the integer part of X, then 

R,(N,X;9)>X + ^^ 

for X > 0. Moreover we have equality whenever the sequence 9 consists of 
equally spaced points. Notice here that X + — ^^)/X > max(l,X) for 
X > 0. 

Part (ii) shows that R{N,X; 9) > X + 0(1) on average with respect to 
X (by virtue of ([5])). Thus the lower bound implicit in Theorem [2] holds, on 
average, for any sequence 9. Moreover, we have 

Ro{N,X^,9) - 1 ^ ^^^^ ^_ ^ ^^^^^ _ 

Hence Theorem [3] shows that for any sequence (6'n)5" one has 

X <i?(X,X;^) <X 
for X > 2, say, providing only that R{N, 1;9) <^ 1. 

The author was introduced to the topic of this article by Jimi Truelsen. 
His input, through a number of interesting conversations, and his helpful 
comments on an earlier draft of this paper, is gratefully acknowledged. 



2 Proof of Theorem 3 



In this section we give the rather easy proof of Theorem [31 For part (i) we 
first show that implies (jl]). We use the fact that 



Ro{N, X-9) = l + - I R{N, t- 9)dt. 



X 



(5) 
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For any fixed positive integer K we have 



k=i 



since -R(A^, t; 6) is non- decreasing with respect to t. We now let tend to 
infinity and apply ([2]), whence 



X V K 



lim sup / R{N, t;e)dt<—y^ 

N-.^ Jo 



X^kX X\K + 1] 



K 2K 



for any fixed positive integer K. Since K is arbitrary it follows that 

lim sup / R{N,t;e)dt < —. 



N^oo Jo 

The corresponding lower bound 



X jr2 

liminf / R{N,t;e)dt > — 







follows similarly from the inequality 

rX V 



/ X kX 
R{N,t-e)dt>-Y,m,^:0). 
k=0 



Thus (ED implies (g]). 

To establish the reverse implication we use a standard Tauberian argu- 
ment. For any A > we have 



R{N,X;9) < / R{N,t;e)dt 

Jx 



= ^^{i?o(iV, X + A; 0) - 1} - ^{RoiN, X- 9) - 1}. 
We let N tend to infinity and apply (jl]) to obtain 

hmsup R{N, X; 9) < ^^(^ + A) - = X + |. 

Since A > was arbitrary we deduce that 

limsupi?(X,X;^) < X 
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The corresponding lower bound is trivial if X = and otherwise follows as 
above, starting with the fact that 



R{N,X;e) > / R{N,t;9)dt 

Jx-A 



= ^{R,iN,X; e)-l}- ^^{i?o(iV,X - A; 6) - 1} 

for any A G (0,X). This establishes part (i) of Theorem El 
For the second part of the theorem we use the fact that 

1 

2, 



where 
However 



XRoiN,X;e)= / Lit.Xydt (6) 
Jo 

L{t, X) = i^{n:\\en-t\\< {2N)-'X}. (7) 

[ L{t,X)dt = [ i^{n:\\9n-t\\<{2N)-^X}dt 
Jo Jo 



N /,e„+(2N)-^X 
'6i„-(2Af)-iX 



n=l 

N 

n=l 

X, 



whence Cauchy's inequality yields 

X^ < [ L{t,Xfdt = XRQ{N,X; 



from which assertion (ii) of Theorem [3] follows, since the terms m = n yield 
RoiN,X;e)>l. 

To handle part (iii) we again use We have 

L{t,X + Y) < L{t+{2N)-^Y,X) + L{t-{2N)-^X,Y) 

^^^ Lit + i2Ny'Y,X) ^ ^^^ L{t-{2N)-^X,Y) 

so that 



by Cauchy's inequality. If we apply © to each term this produces 



(X + Y)MN, X + Y-e)<{X + + ^^^^^^ 

which gives the required inequality. 



3 Proof of Theorem [T] 

In this section we begin by presenting a proof that ([3]) holds for almost all 
real a. This is an immediate consequence of the following result, via the 
Borel-Cantelh Lemma. Of course it suffices to consider a G [0, 1], since the 
set of a for which holds has period 1. 

Lemma 1 There is an explicit sequence of open intervals In, with 

oo 



^^meas(Jn) < cxo, 



n=l 

such that if a & [0, 1] lies in only finitely many of the In then ([3]) holds for 
each fixed X . 

The proof of this result will take up the bulk of this section. It will become 
clear in the course of this work, just what is meant by the word "explicit" in 
the statement above. At the end of this section we shall show how Lemma [1] 
allows us to construct values of a for which holds. 

In the course of the proof we shall use a small parameter rj > 0, which 
we shall take to be 

T] = 1/200 

However we prefer to use the notation rj, which will make it clearer why it 
suffices to use any sufficiently small positive value. The reader will easily 
confirm at each step that rj = 1/200 is indeed suitably small. 

The intervals J„ which we produce will be of three types. We begin by 
including all intervals 

a 1 \ 



q [^2+^"] ' q [q^+^] 
among the /„. Clearly, for these we have 



oo q 



^ meas(/„) ^^^Q ^ ''^J^^ ^ ''<oo. 

n q=l a=0 q=l 
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If a belongs to only finitely many of these intervals we will have 

a 



a 



>q 



-2-7, 



for q > qo{a), say. Then if aj/qj and aj+i/qj+i are successive convergents to 
a, with qj > qo, we will have 



1 



> 



a 



2-ri 



whence qj+i < qj'^^- It follows that if N > qo there will be a convergent a/q 
with 

7V3/2+^ < g < (8) 

and hence 

^2/3-4r,/3 < < q2/3-v/3_ 

Since a/q is a convergent to a we will have 



a 

a 

q 



< 



Let a — a/q + (j), so that |0| < q'"^. We will find it convenient to write a for 
the multiplicative inverse of a modulo q. We begin our analysis of Ra{N, X) 
by observing that 

2R^(N,X) + 1 ^ N-'^{m,n < N : \ \a(m'^ - n^)\\ <X/N}. 

Now, if — = ar(mod q) with |r| < Xq/N — N'^/q, then 

\\a{m^ - n'^)\\ < \\^{m^ - n^)\\ + \(l){m'^ - 71^)1 < \r\q-^ + q-^N^ <XN-\ 

Similarly, if ||Q;(m^ — n^)!] < X/N then — = ar(mod q) with 

|r| < Xq/N + N^q. 

We shall write 

74(A^, g, c) = #{m, n < N : — = c (mod g)}, 

whence 

A^-^ ^ >1(A^, ?,ar) < 2Ra{N,X) + l 

\r\<Xq/N-N'^/q 

< N-^ A{N,q,ar). 

\r\<Xq/N+N^/q 



If we now impose the condition that 

then < Xq/N — N'^/q < Xq/N + N"^ /q < q for large enough N, in view of 
([H]). Moreover, if g f c we have 

A{N,q,c)<^ d{\k\) < N'+'^q-', 

\k\<N^ 
k=c (mod q) 

since the case A; = cannot occur. It follows that 

\r\<Xq/N 

The error term here is 0{N~''^), and 

A{N,q,0) = N + 0{ J2 d{k)) = N + 0{N^+\-^), 

k<N'^,q\k 

whence 

RaiN,X)=N-^ Y A{N,q,ar) + 0{N-''). (9) 

l<r<Xq/N 

We shall see that the expected value of A{N, q, c) is about {N/q)^AQ{q, c), 
where 

Ao{q, c) := #{1 < m, n < q : ~ = c (mod q)}. 

We write 



and 



A(M,g,c) := \A{M,q,c)-{-fAo{q,c) 



A*{q,c) := max A(M,g,c). 

M<q2/3 



With this notation the technical result which is the key to our approach is 
the following. 

Lemma 2 Let q be a positive integer. Write q as a product of prime powers 
in the form q = YlpP'^^^^ ■^^^ 



Qi ■= [[ P 

p=2 or e(p)>l 

Then 



c=l 



with an implied constant which is effectively computable in terms of rj. 
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We will prove this later, in §11 

In view of Lemma [2] we will include among the intervals described in 
Lemma [1] a second category, namely all those 

a 1 a 1 \ 
q q ^^2]y 



for which qi > q^^ . Since 



1/2 



whenever qi > q^'^ we see that for these intervals we have 



oo 1/2 
?1 



^meas(J„)« Y.^'' 



+ + ^1+^ + ^1+^ + • • •} IIU + -j^ + 2^ 1^(1+^ J < 

p>2 ^ e=2 ^ 



From now on we may assume that we have qi < q^"^ for all values of q 
under consideration, so that the estimate in Lemma [2] is of order g3/2+ioj? 
We now define a set B{q) of "bad" values for a by setting 

B{q) ■={0<a<q: ^ A*(g,ar) > g2/3-2^}. 

r<gl/3+2'7 

Then Lemma [2] yields 

q'/'''^#B{q) « E E ^*(^'«^) 

a=0 r<5i/3+»7 



c=l 

by Cauchy's inequality. Thus 

To handle the bad values of a we introduce our third class of intervals /„, 
defined as 

11 



and observe that 

oo 

^meas(/„) < 

n q=l 0<a<q,a€B(q) 

oo 

q=l 

providing that we choose rj < 1/108. Thus for the three classes of intervals 
we have defined we have 

^^meas(/„) < oo 

n 

providing that we choose rj = 1/200, say. 
When a ^ B{q) the estimate ([9]) produces 

Ro.{N,X) = Nq-^ Yl Mq,ar) 

r<Xq/N 

+ 0{N-^ J2 A*{q,ar)) + 0{N-'^) 
= Nq-^ J2 Ao{q,ar) + 0{N-'q^^^-^'^)+0{N-'') 

r<Xq/N 

for g2/3-4r?/3 < jV < g2/3-r?/3^ Under this condition the two error terms are 
both 0{N-i). 

We proceed to investigate 

^/lo(g,ar). 

r<R 

The function Ao^q, r) is multiphcative with respect to q. Thus if qi is defined 
as in Lemma [2] and go '■= we will have A^^q^r) = AQ{qQ,r)AQ{qi,r). 
Since qo is odd we have 

Aoiqo,r) = < go : uv = r (mod go)}, 

whence Ao(go, kr) = y4o(go,r) whenever k is coprime to go. Thus 

91 

^v4o(g, ar) = ^ ^ v4o(go, rf)v4o(gi, go, gi : d, s) 

r<R d\qo s=l 
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where 



U{R,qo,qi : d,s) = 4f {r < R : {qo,r) ^ d, ar = s {mod qi)} 

— ii{e)^{r < R : de\r, ar = s (mod qi)} 

e\qo/d 

- EMe){;^+0(l)}. 

e\qo/d 



It follows that 



say, where 



and 



Y,MQ,ar) = RE + 0{E), 

r<R 



Ql 



E = ^ ^ fi{e)^Ao{qo,d)Ao{qi,s)/deqi 

d\qo e|(}o/rf 
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d\qo e\qo/d s=l 

In a precisely similar way we have 

q 

J2M<l,r)=qE + 0{E), 

r=l 

with the same values of S and E, so that 

R 



r<R ^ r=l 



for R < q. 

Trivially we have 



qi 



^Ao{qi,s) = qf. 



s=l 



Moreover the congruence a;^ = /c(mod go) has ©(gg) solutions a;(mod go), 
uniformly in k, since go is square-free. Thus Ao(go, d) <^ go^^- These bounds 
show that E < go'^^'^gf < g^+^^gi < g^+^''. Finally 

q 

^Ao(g,r) = g^ 
13 



whence 

J]Ao(g,ar) = i?g + 0(gi+^^) 

r<R 

for R < q. We therefore deduce that 

Ra{N, X) = X + OiN-") + 0{Nq"^+^^). 

The final error term will be 0(iV-'') for g2/3-4r,/3 < jy < ^2/3-77/3, 

In conclusion we have shown that holds uniformly for A^"'' < X < N"^, 
providing that a lies in none of the intervals /„. This establishes Lemma [TJ 

To complete the proof of Theorem [1] we describe an algorithm which will 
generate explicit values of a for which ([3]) holds. Suppose we are given a 
closed interval / of positive length, and that we wish to construct a suitable 
a belonging to /. Without loss of generahty we may assume that I C [0, 1], 
since the property (E]) has period 1 in a. Moreover we shall assume that I 
has rational end-points, as we clearly may. Finally we write L for the length 
of I. Now consider the following algorithm. 

We begin by taking rj = 1/200, and we compute an integer such that 

00 

5^ /n < L/2. 

n=N 

We have not specified a numbering for the intervals /„, but it would be easy to 
do so. The contribution from the first two classes of intervals J„ is relatively 
easy to calculate. For the third class one would need to make explicit the 
implied constant in Lemma [2|, but there is no theoretical difficulty in doing 
this. 

Now, for each integer k > N define 

k 

Ffc := / \ U /„. 

n=N 

This is a finite union of closed intervals with rational end points, since the 
intervals J„ also had rational end points. It is important to notice here that 
Ffc cannot be empty, since 

k 

meas(Ffc) > meas(/) — meas(/„,) > L/2 > 0. 

n=N 

We may then compute the set of end points of all the intervals which make 
up Fk, and take to be the smallest such end point. Thus is an explicitly 
computable rational number, with 
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It is clear from the definition that the sets are nested, with 3 
Fn+1 ^ Fiy+2 ■ ■ whence the sequence must be non-decreasing. It follows 
that it converges to a limit, a say. Take any integer j > N. Then, since 
fk^Fk'^ Fj for all k > j, and Fj is closed, it follows that a & Fj. However 
this holds for all j > N, whence 

oo oo 
j=N n=N 

We therefore see that a lies in none of the Ik for k > N, so that holds 
for a, for all X. 

The a G / that we have produced has been "constructed" in the sense 
that we have given a procedure for determining a sequence of rationals which 
converges to a. This completes the proof of Theorem [1] 

4 Proof of Lemma [2] 

We begin this section by considering 

<? 

^:=^A(iV,g,c)^ 

c=l 

for which we prove the following result. 

Lemma 3 Let q be a positive integer, and let qi be defined as in Lemma El 
Then if N < (fl^ we have 

c=l 

with an implied constant which is effectively computable in terms of rj. 
We start by observing that 

c c c 

= Si-2N\-^S2 + N^q~^S^, (10) 
say. Clearly 

Si = X2, X3, 0:4 < : q\xl + - xl - xl}, 
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'S'2 = X3 < A^, X2, X4 < g : q\xl + xj - xj - xj}, 

and 

S3 = i^{xi,X2, X3, XA<q: q\x\ + xl- xl- xl). 

We shall relate Si and 5*2 to ^s, using exponential sums. If we write 
eq{m) := exp{27iim / q) we can use a standard manipulation to show that 

S, = q~^ ^(b; Q)T{bu N, q)T{h2; N, q)T{h; N, q)T{h; N, q), 

b 

where b = (61, 62, 63, ^4) runs over vectors modulo q, and the sums S{h;q) 
and T{b] N, q) are given by 

S{^]q)= 5Z eg(b.x) 

X (mod q) 

and 

T{b;N,q) = ^ e,(-6a;) <min(Ar, 

x<N 

Similarly we find that 

S, = J2 62, 0, 0); q)Tih; N, g)T(62; N, q) 

and 

53 = 5((0,0,0,0);g). 

Since T(0; N,q) = N we see that the terms corresponding to b = cancel in 
fllUp . and it remains to estimate the contribution to Si and S2 arising from 
terms with b 7^ 0. We shall write 

5« = q-^ Yl l)nbi; N, g)T(62; N, q)Tih; N, q)T{h,- N, q) \ 

b 

where the sum is over vectors with \bj\ < g/2, precisely i of which are non- 
zero. Then 

4 

(11) 

i=l 

The sums S{h] q) satisfy a product rule 

-^(b; qiq2) = S{h; qi)S{h; ga), (gi, ^2) = 1 
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and a trivial bound 

\S{h;q)\<q\ 

Moreover wiien q is an odd prime p a standard evaluation shows that 

S{h]p) ^p^+p'^ -p 

when p|b; that S{h;p) = p'^ — p when p\bl + h\ — h\ — h\ but p | b; and that 
S(h]p) — — p if p f 6f + 6| — 6| — 6|. We may therefore decompose q into 
coprimc factors q = qiq2q^qi such that q2, qs, q^ are odd and square-free, with 
q2\h and qslbj + — bl — bl- Moreover we will have 

Sih;q)<^qtql-^\lq,. 

There are 0{q'^) possible factorizations q — qiq2%(lA- Thus 

^««^-4+. qUl-^ylq^S^iq^^q^^q^^q^), (12) 

9l,«2,93,«4 

in which 

SMuq2.q,.qA) = J] |T(6i; iV, g)r(62; iV, g)T(63; iV, g)T(64; iV, | 

b 

4 

< J2ll^HN,q\bj\-'), 

b j=l 

where b runs over vectors in the range \bj\ < q/2, precisely i of which are 
non-zero, and for which 

q2\h, qslbf + bl - bl - bl- 

We shall discuss the case of S^^^ in detail, the other sums being treated 
similarly. For any fixed choice of ± signs we write 

K{Ci,C2,C^,C4;q3) = #{c e : qs\cl±cl±cl±cl \cj\ < Cj, (1 < j < 4)}, 

where we shall assume that 1 < Cj < q for all j. To estimate this we shall 
suppose that Ci < C2 < C3 < C4. If ± c| ± c| ± c| = ^3/0, say, then 
k <^ C^/qs- For each value of k one sees that Ci and C2 determine 0{q^) 
pairs C3, C4, unless cf±C2 = q^k, in which case Ci and C3 determine 0(1) pairs 
C2, C4. Thus there are 0(CiC3g'^') possibilities for each value of k, so that 

KiCi, C2, Cs, C4; qs) « (1 + Cl/qs)CiCsq^. 
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For an alternative estimate we observe that the congruence c1 = ?T,(mod q^) 
has 0(^3) solutions modulo q^, whence 

K{C^, C2, C,, qs) « (1 + CJq,)C^C2Csq''. 

To put these bounds into a more convenient form we write C := C1C2C3C4 
and observe that by combining our estimates we have 

K{Ci, C2, C3, C4; ga) « {C1C3 + q^'C + mmiq^'CC^Cl , CiCaCa)}?''. 
Since Ci < C2 < C3 < C4 we have C1C3 < VC and 

mm{q^'CiCsCl,CiC2C3) < {q^'CiCsClY/\C,C2CsY/^ 

— ^3 ^v^2 ^'y^^ 

where we have written Cq = minCj. It follows that 

i^(Ci,C2,C3,C4;g3) « {C^'^ + q^^C + q^''^CC,"^}q^ 

« C{C^' + q^' + q-'''C,''']q\ 

We can now bound 5*4(^1, Q's, Q'4)- We write hj = q2Cj and decompose 
the range for each Cj into intervals either of the shape \cj\ < Cj = gg^^A^"^, 
or of the form Cj/2 < \cj\ < Cj with Cj > qq^'^N'^. Ranges with Cj < 1 
will not arise, since all hj are non-zero for S^'^\ There will be <C log q <^ q^ 
sets of ranges in total, on each of which we will have 

4 

J]min(iV,g|6,|-i)«gV'C"'- 
i=i 

Moreover, since Cq > qq2^N^^ > q^^^q2^, we find that 

K{CuC2,C3,C,;qs) « C{q-'/'ql + q^' + q,'/\~'/'ql/'}q^ 

« CqtV\^- 

It now follows that 



54(gi, 52,53,94) « q''-''q2'C-\Cqt^\-'^\^ 
whence (fT2l) yields 



<l'^''<l2V 



« q'/'-^'^ql 
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Similar arguments show that 

#{c G I? : q^\c{ ±c\± c^, |c,| < C„ (1 < J < 3)} 

whence 

5(3) iVg2/3+4.g3 <^ ^4/3+4,^3. 

that 

#{c G Z2 : q^\c\ ± c^, |c,| < C„ (1 < J < 2)} « + g3-^C2)g3^ 

« (ClC2)g2'^'g3^'^\^ 

whence 

5(2) « iV2g-l/3+4r,^3 ^1+4,^3. 

and that 

#{cGZ:g3|c^ |c| < C} < C/g3 

whence 

5(1) « iV3g-l+4''g3 ^l+4,^3_ 

In view of ( fTTj) these estimates suffice for the proof of Lemma [3l 

We proceed to deduce Lemma [2] from Lemma [31 If < M < then 

A{N + M, g, c) - A{N, q, c) 

■C #{m G (A^, + M], X < + M : - = ±c (mod g)}, 

whence 

q 

Y,{A{N + M,q,c)-A{N,q,c)Y 

c=l 

< t; G (A^, A^ + M], X, 2/ < A^ + M : q\u^ - ± {v^ - y^)}. 

If — x^ ± (f 2 — y^) = kq <^ N"^ then each set of values u, v, k determines 
0{q^) pairs x,|/, except when the ib sign is negative and — f ^ = kq. 
Thus there are 0{M'^N'^q~^^^) solutions u,v,x,y with — ^ kq. When 
— f 2 = kq we have k <^ MNq~^. Thus there are 0{MNq^^^^) pairs m,!" 
corresponding to non-zero values of /c, and M pairs for = 0. To each such 
pair u, V with u'^ — v'^ = kq there correspond 0{N) pairs m = n. We therefore 
obtain the bound 

J]{A(Ar + M, g, c) - y4(Ar, g, c)}^ < M'^N'^q-^+i + MN'^q-^+'^ + MAT. 

c=l 
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In particular, taking M = N = q, we have 



c=l 

It follows that 

' 0<H<M 



c=l 

9 



« 5Z max {A(iV + i7,g,c)-A(iV,g,c)}^ 



' 0<H<M 
c=l ~ ~ 



+ max '-^Aoiq, 

0<H<M q'^ ^ 

c=l 



^{A(iV + M,g,c)-A(iV,g,c)}^ 



c=l 



c=l 

< ^''{M^A^^^-^ + MAT} 



if AT < g2/3 
To handle 



A*(g, c) = max A{U,q,c) 



we cover the available range for U with 0{q^^^M ^) sub-intervals Nj < U < 
Nj + M, whence Lemma |3] yields 



,0^ 



c=l 



1 

« 5^ J2^A{N„ q, cf + max \A{N, + if, g, c) - A(Ar,, g, c) 1^} 



0<H<M 

Ni c=l - - 



<^ g2/3^-1^4,|^4/3^3 ^ ^1/3^2 ^ ^2/3^|^ 

The choice M = g^/^ then results in the estimate required for Lemma [21 
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5 Proof of Theorem [2 

On writing n — m = u, n + m = fwe find that 

Ra{N,X) = N-^#{m <n<N : \\a{m'^ - n^)\\ < XN''^} 
= N-^ = u (mod 2) : u < v < 2N - u, \\auv\ \ < XN'^}. 

u<N 

When u is even we write 2u in place of u and put f = 2x to find that the 
corresponding contribution is 

AT-i #{a;GN: u < x < N - u, \\Aaux\\ <XN-^} 

u<N/2 

= iV"^ J] {i?(Ar-u,4au,X/Ar) -i?(u,4au,X/Ar)}, (13) 

u<N/2 

say, where 

R{M,P,S) := #{x G N : x < M, \\px\\ < 5}. 

To handle odd values u we count integers v = 1 (mod 2) by first considering 
the contribution from all v, and then subtracting the contribution from even 
V. This leads to a total 

Y {R{2N -u,au,X/N) - R{u,au,X/N)} 

u<N, 2\u 

-N-^ J2 {R{N -u/2,2au,X/N) - R{u/2,2au,X/N)}. 

u<N, 2\u 

We proceed to estimate i?(M, /3, 5). Let u, v G be the vectors 

Then for 5 G (0, 1) we have 

1 + 2i?(M, /3, 5) = #{(x,i/) G : |x| < M, < 5} 



= #{(x,i/) G : xu + yv G [-VM5, VM5]2}. 
The vectors u, v generate a lattice of determinant 1 . Hence 

#{(x, G Z2 : xu + G [ -S, S ]''} = {2Sf + 0(8/ \^) + 0(1), (14) 
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where Ai is the first successive minimum of the lattice, that is to say the 
length of the shortest non-zero vector in the lattice. In our case we find that 

1 + 2R{M, (3, 6) = 4MS + 0{M^/^6^/^X^^) + 0(1). 

We have 6 different pairs (M, f3) = {N — u, 4au), . . . ,{u/2, 2au) to consider, 
each with a corresponding value for Ai. We write Aq for the smallest of these 
6 values, and split the available range for M^/^(5^/^Ao ^ into dyadic intervals 

AE < M^/^^^/^Ao ^ < 8E. 

For values u for which E < X^^'^ the total contribution of the error terms to 
Ra{N,X) is clearly The choice of X^/^ as the point at which we 

split the range for E is not optimal, but is adequate for our purposes. 

In the remaining case E > X^/^ there will be coprime integers x, y for 
which \x\ < M/{AE) and \l3x - y\ < 5/{AE). It follows that we will have 

R{N -u,Aau,X/N) = 2X{1 - u/N) + 0{E), 

R{u,Aau,X/N) = 2Xu/N + 0{E), 

R{2N -u,au,X/N) = 2X{2 - u/N) + 0{E), 

R{u,au,X/N) = 2Xu/N + 0{E) 

R{N -u/2,2au,X/N) = X{2 - u/N) + 0{E), 

R{u/2,2au,X/N) = Xu/N + 0{E), 

and that there is a coprime pair x, y satisfying one of 

\x\ < ——, \Aaux — y\ < 



or 



or 



AE' ' - ANE 



N . . X 



I I - 2£;' ' ^' - ANE 

, , iV , , X 

\x\ < —-=;, \2aux — y\ < 



AE' ' - ANE 

respectively. To simplify matters we replace xhy x' = Ax in the first case and 
by x' = 2x in the third, and then remove a factor {x', = 2 or 4 if necessary. 
We deduce in each case that there is a coprime pair with \x\ < N/E and 
\aux — y\ < X/ {NE). Since they are coprime, x and y cannot both vanish. 
Indeed, since X < log N we will have X/ (NE) < 1 whence it is clear that x 
cannot vanish. It follows that the total contribution of the error terms 0{E) 
to Ra{N,X), arising from an individual value of £^ > X^/^, is 

E N X 

< — x,y)eN^ xZ: u<N,x<—,{x,y) = 1, \aux - y\ < -^}. 



22 



We now calculate the contribution from the main term of R{M, (3, 6). For 
(fT^ this is 

1 AX X X 

- Y: {2^(1 - ^/N) - 2XU/N} = ^ E (^/2 - = y + 0{j^). 

u<N/2 u<N/2 

Similarly the odd values of u contribute X/2 + 0{X/N). We therefore con- 
clude as follows. 

Lemma 4 If 1 < X < log N then 

RaiN,X)=X + OiX^/^) + I y EV*{N,^; ^ 



X^/^<E=2''<N 

where 

V*{A,B;A) : = 

X, ?/) e X Z : u < A, X < B, (x, y) = 1, \aux - y\ < A}. 

It is already clear here that our approach cannot provide an asymptotic 
evaluation for Ra{N,X) unless X oo. The error term 0(1) in (IT^ will 
produce at least a corresponding error 0(1) for Ra{N,X). Any sharper 
estimate would appear to require information on the way the shape of our 
lattice varies with the parameter u. 

From now on we shall focus on the second error term above. We write 
(m, y) = f and suppose that / lies in a dyadic range F < f < 2F. Given such 
an f, a u = fuo and y = fyo then uq < N/F and \auoX — yo\ < X/{NEF). 
Moreover each pair uq, yo can correspond to at most F pairs u,y, since we 
are assuming that F < f < 2F. Our error term is therefore 

«iV-i y y EFVi-,-;^) (15) 

^e'f'nef' ^ ' 

Xi/2<£;=2'=<Af l<F=2''<Ar 

where 

\/(A5;A) := 

#{(a, 6, z) e X Z : a<A,h< B, {ab, z) = 1, \aab - z\ < A}. 
Our strategy for tackling V{A, B; A) is based on the following lemma. 

Lemma 5 Let E be an ellipse centred at the origin, of area A{E). Then the 
number of coprime integer pairs [x, y) E E is 0(1 + A{E)). 
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This easy result may be found in the author's work [21 Lemma 2], for example. 
If we fix a, say, then the ellipse 

{{b, z)eR^ : + A->a6 - < 2} 

has area -C BA, and we deduce that V{A, B; A) ^ ABA-\-A. By symmetry 
we then have the bound 

V{A,B;A) <^ABA + mm{A,B). (16) 

In our application the contribution from the term ABA is usually satisfac- 
tory, but the effect of the second term is likely to be too large unless A = N/ E 
and B = N/ F have very different sizes. To circumvent this difficulty we shall 
use a delicate arithmetic trick, which is the key to our attack on Theorem [21 
It will be convenient to assume that A < B, a.s we may, by symmetry. 
We take parameters -Pi > -Po > 1 and consider prime factors p of ab in the 
range Pq < p < Pi- Thus we will need to consider separately 

Vi ■= #{(a, 6, ;2) G X Z : a < A, b < B, {ab, zU) = 1, \aab - z\ < A}, 

where 

n:= Hp 

Po<P<Pi 

Here Lemma [S] shows that 

Vi < J2 #{{b,z) eZ^ : b< B, {b,z) = 1, \aab- z\ < A} 

a<A,{a,n)=l 

< {BA + l)#{a <A: {a,U) = 1}. 

The number of available integers a may be estimated using a standard sieve 
bound. According to Theorem 2.2 of Halberstam and Richert [T], for exam- 
ple, one has 

#{a < A : (a,n) = 1} « A J] - Vp) « 

Po<p<Pi ^ 

providing that Pi < A. This yields the following lemma. 
Lemma 6 If 1 < Pq < Pi < A we have 

V,«ABA + A^:^. 

log Pi 
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If ab does have a prime factor p in the range Pq < p < Pi we may 
choose the smallest such prime p, and classify the corresponding triples a, b, z 
according to the dyadic range P = 2^ < p < 2P in which p lies. We write 
V2{P) for the corresponding contribution to V{A, B; A) and write 

V, := Yl 

Po<P=2''<Pi 

If p\a we set a' = a/p and b' = bp, while if p f a we will have p\b, and we 
set a' = ap and b' = b/p. It follows that \aa'b' — z\ < A, and that either 
a' < A/P, b' < 2BP or a' < 2AP, b' < B/P. Moreover, if we are given 
a triple a',b',z counted by V{A/ P,2BP\ /\) then it determines the prime 
p, which will be the smallest prime p > Pq dividing a'b' = ab. Knowing p 
one may then find the pair a, b which must either be a = a'p, b = b' /p or 
a = a' /p, b = b'p. It follows that each triple a', b', z counted by 

V{A/P, 2BP; A) + V{2AP,B/P; A) 

arises from at most 2 triples a, b, z counted by V2{P). We may now use (fT6l) 
to deduce that 

V2{P) < ABA + mm{A/P, BP) + mm{AP,B/P) < ABA + B/P, 
from which we obtain the following lemma. 
Lemma 7 If 1 < Pq < Pi we have 

V2<^ABA(\ogPi) + Po^B. 



We now use Lemmas [6] and [7] to estimate the contribution to (fT5l) from 
terms with EF < (logA^)^/^. We choose 

Pq = E^F\ Pi = Qexp{{EFf^), 

so that 1 < Pq < Pi < A for large enough N. The terms ABA and 
ABA{logPi) in Lemmas [6] and [7| then contribute 



« y y EF^-f^{EFf' 

^ ^ EFNEF^ ' 

X^/^<E=2''<N 1<F=2'^<N 

^ Yl Yl jefw 

« x^/^ 
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Moreover the error term A{log Pq) / {log Pi) in Lemma [U] produces 

since min(A^/E, A^/F) < N{EF)-'^'^. Finally, the error term Pq occurring 
in Lemma [7] produces 

«^"' E E i?^max(|,|)i?-^F-^«L 

X^/'^<E=2^<N l<F=2h<N 

Thus those terms with EF < {\ogN)^^^ make a satisfactory contribution in 
Theorem [21 

Up to this point we have made no use of the Diophantine approximation 
properties of a, but it is time to bring these into play. In order to clarify 
the rationale behind our choice of the various exponents which will occur, we 
introduce constants /3, 7 G (0, 1) on which we will impose certain constraints 
as the argument progresses, and which will eventually be specified in 
To begin with we assume that a satisfies 



a 

a 



1 

> 



for every fraction a/q. In particular, if V{N/ E, N/ F; X/{NEF)) counts 
x,y,z, so that \a — z/xy\ < X/{NEFxy), we deduce that 

{xyY^^ > NEF/{kX). 

Hence V{N/ E, N/ F; X/{NEF)) = unless 

{N^/EFf^^ > NEF/{kX). 

We therefore assume from now on that 

{logNf/^ <EF< (^X)l/(2+/3)^{l+2/3)/(2+/3)_ (^7) 

Let 

Q = 



EF 

and apply Dirichlet's Approximation Theorem to obtain coprime integers a, q 
with 

a 1 ^ 
a-- <—, l<q<Q. 

q qQ 
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It follows of course that 

(/?-ig)i/(i+^) <q<Q. (19) 
Now if \\axy\ \ < X/{NEF) with x < N/E and y < N/F then 

\\axy/q\ \ < X/{NEF) + N'^/{EFqQ). 
whence axy = r(mod q) for some integer r with 

|r| < qX/{NEF) + N'^/{EFQ). 

It follows that 

\r\<qX/(NEF)+N2/{EFQ) n<N^/{EF) 

an=r (mod q) 

where d{n) is the divisor function. The reader should observe that there is a 
loss at this point, in replacing 

y : X < N/E, y < N/F, xy = n} 

by d{n). This loss is of order log A^, and is only acceptable since we are now 
in the case in which EE is of larger order than log A^. 

For the case in which [s, g) = 1 it was shown by Linnik and Vinogradov 
[3] that 

y ci(m) log M, 

m<M ^ 
m=s (mod q) 

providing that q < M^^'^ for some constant 7 > 0. In general, if {s, q) = h 
say, one may write s = hs' , q = hq' and m = hm', so that q' < {M/hY^'^ 
and (i(m) < d{h)d{m'). Then 

J2 d{m) < d{h) J2 ^K) 

m<M m'<M/h 
m=s (mod q) rn'=s' (mod q') 



Q 

<^ d{h)q'^M\ogM. 



«, d{hf-^'-^\ogM 
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If we sum for Isl < S* we find that 



J2diis,q)) < \s\ < S, k\s} 

\s\<S k\q 

< ^d{k){l + S/k) 

k\q 

« d\q) + S\{{l + 2/p + 0{p-^)) 



< d\q) + S 



g2 



2 



< d'^{q) + Si\og\ogqf. 

It follows from ([20]) that 

^E' F' NEF' 

„v T\j2 t\t2 

« KW + (^ + ^)(loglog,n^log]V. (21) 

We begin by examining the first term d'^{q)N'^{qEF)~^ log A^. In view of 
ffT^ and (dHD this is at most 

at2 / Ar2 \ {/3+7)/(l+/9) 

/s:i/(i+^)d2(g)(logiV)4I-g-i/(i+^) = K^/(i+^)d2(g)(logiV) ^ ^ 



■■) 



EF^ v^/v o . y^^^ 

When we multiply by N~^EF and sum over dyadic ranges subject to ffT7|) 
we see that the contribution to ( fTSl) is 

with 

01 = 1 - (1 + /5)-M2(/? + 7) + (1 - 7)\^}. 

Turning to the second term on the right of (1211) . we see that the overall 
contribution to the error terms in Theorem [2] is 

X (log log g) 2 (log A^) {EF)-^ <^X{\og\og Nf {log N)'^/\ 

-EF>{logAf)5/4 

which is 0(X^/®). This is satisfactory for the theorem. 
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Finally, the third term on the right of (^T^ is 

« (loglogiV)2(logiV) (^^y «:i/(i+/3)g-(2+/3)/(i+/3) 

2 N (27+/37+/3)/(l+/3) 



(loglogiV)2(logiV) i^—j 



We impose the condition that /5,7 < 1/3, which ensures that the exponent 
(27 + jSj + l3)/{l+ (3) is less than 1. Now, when we multiply by N~^EF and 
sum over dyadic ranges subject to ffT7|) . we get an overall contribution 

« {kX) (1-27-^«7)/(1+/3){2+/3) (log log iV)2 (log N)'N-''>' , 

with 



If 7 were equal to zero we would have 

1-3(3 -(3"^ 



(l + ^)(2 + /3)- 

We will first choose (3 so as to make this value positive, and then select a 
sufficiently small 7 so that 0i, 02 > 0. With this in mind we specify 

/3=l/4, 7 = 1/40, (22) 

from which the assertion of Theorem [2] follows. 
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